Bell inequalities for arbitrarily high dimensional systems 
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We develop a novel approach to Bell inequalities based on a constraint that the correlations 
exhibited by local variable theories must satisfy. This is used to construct a family of Bell inequalities 
for bipartite quantum systems of arbitrarily high dimensionality which are strongly resistant to noise. 
In particular our work gives an analytic description of previous numerical results, and generalizes 
them to arbitrarily high dimensionality. 



One of the most remarkable aspect of quantum me- 
chanics is its predicted correlations. Indeed, the cor- 
relations between outcomes of measurements performed 
on systems composed of several parts in an entangled 
state have no classical analog. The most striking aspect 
of this characteristic feature of quantum physics is re- 
vealed when the parts are spatially separated: no clas- 
sical theory based on local variables can reproduce the 
quantum correlations. Historically, this became known 
as the EPR paradox and was formulated in terms of mea- 
surable quantities by Bell and by Clauser, Horn, Shi- 
mony and Holt Q as the nowadays famous inequalities. 
Other aspects of quantum correlation were analyzed in 
the form of paradoxes, like, e.g. Schrodinger's cat and 
the measurement problem. In recent years, these para- 
doxical aspects have been overthrown by a more effective 
approach: let 's exploit "quantum strangeness " to perform 
tasks that are classically impossible has become the new 
leitmotiv! From this "conceptual revolution" , the field of 
quantum information emerged. Old words became fash- 
ionable, like "entanglement" . Old questions were revis- 
ited, like the classifications of quantum correlation. 

The variety of known partial results, in particular 
about entanglement measures, makes it today obvious 
that there is no one-parameter classification of entan- 
glement. This letter concerns classifications related to 
what is called quantum non-locality, i.e. the impossi- 
bility to reproduce quantum correlations with theories 
based on local variables (often called local realistic theo- 
ries). Specifically we develop a powerful new approach to 
Bell inequalities which we then use to write several fam- 
ilies of Bell inequalities for higher dimensional systems. 

Local variable theories cannot exhibit arbitrary corre- 
lations. Rather the conditions these correlations must 
obey can always be written as inequalities (the Bell in- 
equalities) which the joint probabilities of outcomes must 
satisfy. Our approach to Bell inequalities is based on a 
logical constraint the correlations must satisfy in the case 
of local variable theories. In order to introduce this con- 
straint, let us suppose that one of the parties, Alice, can 
carry out two possible measurements, Ai or A2, and that 
the other party. Bob, can carry out two possible mea- 



surements, Bi or B2. Each measurement may have d 
possible outcomes: Ai, A2, Bi, B2 = 0, . . . ,d ~ 1. With- 
out loss of generality a local variable theory can be de- 
scribed by probabilities Cjkim {j,k,l,m = 0, . . . ,d— 1) 
that Alice's local variable (jk) specifies that measure- 
ment Ai gives outcome j and that measurement A2 gives 
outcome k; and that Bob's local variable {Im) specifies 
that measurement Bi gives outcome I and that measure- 
ment B2 gives outcome m. (In this formulation Alice 
and Bob's strategy is deterministic since it is completely 
determined by the value of their variables jk and Im. 
Any non deterministic local theory can be rephrased in 
the above way by incorporating the local randomness in 
the probabilities Cjkirm see for instance |^). Since they 
are probabilities the Cjkim are positive {cjkim > 0) and 
sum to one (J^jkim'^jkirn — !)• The joint probabilities 
take the form P{Ai = j, Bi ^ I) = J2km ^jkim, and sim- 
ilarly for P{Ai = j,B2 = m), P{A2 = fc, Si = /) and 
P(^2 = k,B2^ m). 

Let us consider a particular choice of local variables 
jklm (this choice occurs with probability Cjkim)- Since 
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We see that the difference, r', between Ai and Bi can 
be freely chosen by choosing j and I. Similarly the dif- 
ference, s', between Bi and A2 and the difference, t', 
between A2 and B2 can be freely chosen. But then the 
difference u' between B2 and Ai is constrained since we 
necessarily have 

r' + s' + t' + u' = . (2) 

Thus in a local variable theory the relation between three 
pairs of operators can be freely chosen, but then the last 
relation is constrained. 

This constraint plays a central role in our Bell inequal- 
ities. Indeed they are written in such a way that their 
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maximum value can be attained only if this constraint is 
frustrated. The simplest such Bell expression is 



+P{A2 = B2) + P{B2^Ai) 



(3) 



where we have introduced the probability P{Aa — Bb+k) 
that the measurements Aa and Bi, have outcomes that 
differ, modulo d, by k: 



d-l 



P(Aa =Bb + k) = J2 Pi^a = j, Bb^j + k mod d) 
j=o 



(4) 



Because the difference between Aa and Bb is evaluated 
modulo d, all the outcomes of Aa and Bb are treated on 
an equal footing. As we see in eq. (|^) this symmetriza- 
tion is the key to reducing Bell inequalities to the logical 
constraint that is imposed by local variable theories. In- 
deed because of the constraint eq. (||) any choice of local 
variables jklm can satisfy only three of the relations ap- 
pearing in eq. (||), eg. Ai = Bi, Bi — A2 etc. . . . 
Hence /(local realism) < 3. On the other hand non-local 
correlations can attain 1 = 4 since they can satisfy all 4 
relations. 

In the case of two dimensional systems the inequality 
/(local variable) < 3 is equivalent to the CHSH inequal- 
ity But the power of our reformulation is already 
apparent since this inequality generalizes the CHSH in- 
equality to arbitrarily large dimensions. In fact the above 
formulation of the constraint imposed by local realistic 
theories allows one to write in a unified way all previously 
known Bell inequalities It can also serve to write 
completely new Bell inequalities and this is the subject 
of the present article. Specifically we have generalised in 
a non trivial way (see eqs. (||) and (||) below) the Bell 
expression to d dimensional systems (for any d> 2). 

One of the interests of these new Bell expressions is 
that they are highly resistant to noise. Indeed Bell in- 
equalities are sensitive to the presence of noise and above 
a certain amount of noise the Bell inequalities will cease 
to be violated by a quantum system. However it has been 
shown by numerical optimization that using higher 
dimensional systems can increase the resistance to noise. 
The measurements that are carried out on the quantum 
system in order to obtain an increased violation have 
been described analytically in Q. And an analytical 
proof of the greater robustness of quantum systems of 
dimension 3 was given in |^]. One of the interests of 
our new Bell inequalities is that when we apply them 
to the quantum state and measurement described in |^] 
for those dimensions {d < 16) for which a numerical op- 
timisation was carried out in Q, we obtain the same 
resistance to noise as in 

The first generalisation of the Bell expression eq. (|^) 

is 

/3 = + [PiAi = Bi) + P{Bi =A2 + 1) 



+P{A2 = B2) + P{B2 = Ai)] 
[P{Ai = Bi - I) + P{Bi = A2) 
+P{A2 = B2-l) + P{B2 = Ai 



1)] • (5) 



The maximum value of Ij, for non-local theories is 4 
since a non local theory could satisfy all 4 relations that 
have a -I- sign in (^). On the other hand for a local vari- 
able theory Iz < 2. This should be compared to the 
constraint /(local variable) < 3 for the expression (||). 
The origin of this difference is the — signs in (|). Indeed 
we have seen when analyzing (^) that only three of the 
relations with a -I- sign can be satisfied by local realistic 
theories. But if 3 relations with + are satisfied in (H), 
then necessarily one relation with — is also satisfied giv- 
ing a total of /a = 2. Alternatively one can satisfy 2 
relations with + and two relations with weight zero (if 
the dimension is larger than 2), once more giving a total 
of h = 2. 

For d = 2 the inequality /s (local variable) < 2 is equiv- 
alent to the inequality /(local variable) < 3 and therefore 
to the CHSH inequality. But for d > 3 the inequality 
based on I3 is not equivalent to that based on /. For the 
quantum measurement described below (when d > 3) the 
inequality based on I3 (and its generalisations Id given 
below) is more robust than that based on /. 

The Bell expression I3 can be further generalised 
when the dimensionality is greater than 3 by adding 
extra terms. The extra terms in Id do not change 
the maximum value attainable by local variable theories 
(/™°^ (local variable) = 2), nor do they change the max- 
imum value attainable by completely non local theories 
^jmax _ ^-j^ However these extra terms allow a bet- 
ter exploitation of the correlations exhibited by quantum 
systems. 

These new Bell expressions have the form: 
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( + [P{Ai =Bi+k) + P{Bi = A2 + k+l) 

+P{A2 = B2 + k) + P[B2 = Ai + k)] 
- [P{Ai =Bi-k-l) + P{Bi =A2-k) 

+P{A2 = B2-k-l)+ P{B2 =Ai-k- 1)]) . (6) 

As mentioned above the maximum value of Id is 4. This 
follows immediatly from the fact that the maximum 
weight of the terms in (^ is -1-1. And the maximum 
value of Id for local variable theories is 2. We now prove 
this last result. 

The proof consists of enumerating all the possible rela- 
tions between Ai, Bi, A2, B2 allowed by the constraints 
(^). This is most easily done by first changing notation. 
We do not use the coefficients r',s',t',u' defined in (|^), 
but use new coefficients r, s, t, u defined by the relation 



2 



Ai= Bi+r , Bi= A2 + S + I , 
= S2 + t , B2= Ai+u , 

which obey the constraint 

r + s + t + u+1 = mod d . 



(7) 



(8) 



Furthermore we restrict (without loss of generahty) 
r, s, t, u to he in the interval 



[d/2] <r,s,t,u< [(d-l)/2] 



(9) 



With this notation the value of the Bell inequality for a 
given choice of r, s,t,u is 

/rf(r, s, t, u) = f{r) + f{s) + fit) + f{u) (10) 

where / is given by 



(11) 



We now consider different cases according to the signs of 
r, s, t, u. 

1. r,s,t,u are all positive. Then (||) and (|) imply 
that r + s + t + M = d— 1. Inserting into (1^ and 
using (|ll|) one finds Id = 2. 

2. Three of the numbers r,s,t,u are positive, one is 
strictly negative. Then (||) and (^) imply that ei- 
ther r + s + t + u — d— 1 oir + s + t + u = —1. 



Inserting into ( |10[ ) and using (11) one finds either 
Id = -2/{d-l) OT Id = 2. 

3. Two of the numbers r,s,t,u are positive, two are 
strictly negative. Then (||) and (||) imply that 
r + s + t + u — —1. Inserting into ([To|) and us- 
ing (|l|) one finds Id = -2/{d - 1). 

4. One of the numbers r,s,t,u is positive, three are 
strictly negative. Then (||) and (^) imply that ei- 
ther r + s + t + u = —1 or r + s + t + u = —d — 1 . 
Inserting into (|lO[) and using ( pT] ) one finds either 
/d = -2(d+l)/p-l) or/d = -2/(d-l). 

5. The numbers r, s, u are all strictly negative. Then 
(1) and (1) imply that r + s + t + u = -d - 1. 
Inserting into (10|) and using ( |ll|) one finds = 
-2(d+l)/(d- 1). 

(Note that for small dimensions d not all the possibilities 
enumerated above can occur. For instance for d = 2, the 
only possible values are Id — ±2.) Thus for all possible 
choices of r, s, t, u, /^(local realism) < 2. This concludes 
the proof. 

Let us now consider the maximum value that can be 
attained for the Bell expressions Id for quantum measure- 
ments on an entangled quantum state. We have carried 
out a numerical search for the optimal measurements. It 
turns out that the best measurements that we have found 



numerically give the same value as the measurements de- 
scribed in . We do not have a proof that these mea- 
surements are optimal, but our numerical work and the 
numerical work that inspired 1^ suggests that this is the 
case. 

We therefore first recall the state state and the mea- 
surement described in |^ . The quantum state is the max- 
imally entangled state of two d-dimensional systems 
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The measurements is carried out in 3 steps. First Al- 
ice and Bob give each of the states \j) a variable phase, 
gi^aU) for Alice and e^'^''^^^ for Bob, which depends on 
the measurement they want to carry out. The state thus 
becomes 
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\j)a ® \3)b 



(13) 



where 0i(j) = f a^, 02(j) - f «2j, ^y{3) = f 
and 1/72 (j) = ^/^2j with ai = 0, a2 = 1/2, /3i = 1/4 and 
/32 — —1/4. The second step consists of each party car- 
rying out a discrete Fourier transform to bring the state 
to the form 
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The final step is for Alice to measure the k basis and Bob 
to measure the I basis. Thus the joint probabilities are 



PQM{Aa ^k,Bb^ I) 
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^ly^exp 
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(3b) 



1 sin^ [7r(fc - I + aa+ Pb)] 
d^ sin^[7r(fc ~ I + Ua + l3b)/d] 



2d^ sin^[7r(fc ~ I + Ua + Pb)/d] 



(15) 



where in the last line we have used the values of aa and 
Pb given above. 

Equation (|l^) shows that these joint probabilities have 
several symmetries. First of all we have the relation 

PQM{Aa =k,Bb = l)= PQM{Aa = k + c, Bb = I + c) 

for all integers c. This symmetry property justifies us 
considering, as in (Q), only the probabilities that Aa and 
Bb differ by a given constant integer c: 

PQAliAa = Bb + c)=Y^ PQAliAa = j + C, Bb ^ j) 



= dPQMiAa=C,Bb = 0) 



(16) 
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Furthermore we have the relation 



Pqm{Bi=A2 + c+1) 



Pqm{A2 =B2+c) = Pqm{B2 =Ai+c) 



(17) 



Using eqs. (|15| to 17) we can rank these probabihties by 
decreasing order. Let us denote 

Qc = Pqm{Ai = Bi + c) = 1/ (2^3 sin2[7r(c + l/4)/d]) . 

(18) 



Then we have 

qa > q-i > qi > q-2 > 92 > 



> Q-ld/2] (> q[d/2]) (19) 



where [x] denotes the integer part of x and the last term 
between parenthesis occurs only for odd dimension d. 
This suggests that the quantum probabilities violate the 
constraints imposed by local variable theories. Indeed the 
probabilities in ( pT| ) are maximized by taking c = 0, but 
then the 4 relations that appear in (|l7| ) are incompatible 
with local realism. In fact replacing the above probabil- 
ities in the expression (^) yields a value Iqm — 4(iqo > 3 
for all dimensions d. 

However a stronger violation is obtained if instead of 
using the Bell expression J, one uses the Bell expressions 
Id- In fact for a d dimensional quantum systems, one can 
use all the Bell expressions Ik for k < d, but the strongest 
violation is obtained by using the Bell expression Id- This 
value, denoted Id{QM), is given by 
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(20) 



For instance we find 

IsiQM) = 4/ (-9 + eVtj ~ 2.87293 , 

h{QM) = ^ (^V2 + y/lO~V?j ~ 2.89624 , 

1 



lim IdiQM) = 



1 



(fc + 1/4)2 (fc + 3/4)2 



= 32 Catalan/vr^ ~ 2.6981 

where Catalan ~ 0.9159 is Catalan's constant. 

In the presence of uncolored noise the quantum state 
becomes 



X 1 



(21) 



where p is the probability that the state is unaffected by 
noise. The value of the Bell inequality for the state p is 

Id{p) ^ pIdiQM) (22) 

Hence the Bell inequality Id is certainly violated if 
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(23) 



(If there is a quantum measurement giving a value of Id 
greater than that given by eq. (|20|), then of course the 
Bell inequality would be violated with even more noise. 
This remark applies to the various p™™ below). 

As a function of d one finds that p™™ is a decreasing 
function of d. For instance: 
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0.69615 

- 72) ~ 0.69055 



lim ""V (16 Catalan) ~ 0.67344 

d — >oo 

For 0? = 3 this reproduces the analytical result of . And 
combining eqs. (^ ) and (23) reproduces the numerical 
results of |^ for all dimensions (2 < d < 16) for which a 
numerical optimization was carried out. 

In summary our reformulation of Bell inequalities in 
terms of a logical constraint local variable theories must 
satisfy has provided the basis for constructing a large 
family of Bell inequalities for systems of large dimension. 
The numerical work of |,§ and a numerical search of 
our own suggest that these Bell inequalities are optimal 
in the same sense that the CHSH inequality is optimal 
for 2 dimensional systems, namely both the resistance 
to noise and the amount by which the inequality are vi- 
olated are maximal. For this reason we hope that the 
Bell inequalities presented here will have as much inter- 
est for physicists studying entanglement of systems of 
large dimensionality as the CHSH inequalities have had 
for bidimensional systems. 

Note: while completing this paper we learned of a Bell 
inequality for qutrits |8j that exhibits the same resistance 
to noise as that obtained in 1^-0]. 
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